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ABSTRACT 

The role of streamwise-spanwise vorticity interactions that occur in turbulent shear flows on flame/vortex 
interactions is examined by means of asymptotic analysis and numerical simulation in the limit of small 
Mach number. An idealized model is employed to describe the interaction process. The model consists of a 
one-step, irreversible Arrhenius reaction between initially unmixed species occupying adjacent half-planes 
which are then allowed to mix and react in the presence of a streamwise vortex embedded in a shear flow. 
It is found that the interaction of the streamwise vortex with shear gives rise to small-scale velocity oscilla- 
tions which increase in magnitude with shear strength. These oscillations give rise to regions of strong tem- 
perature gradients via viscous heating, which can lead to multiple ignition points and substantially decrease 
ignition times. The evolution in time of the temperature and mass-fraction fields is followed, and emphasis 
is placed on the ignition time and structure as a function of vortex and shear strength. 
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1. INTRODUCTION 


The presence of streamwise vortices is now very well established in wall-bounded shear flows as well 
as free shear flows either in transitional state or fully developed turbulent state. The investigations of 
Klebanoff, Tidstrom and Sargent (1962), Wray and Hussaini (1984), and Zang and Hussaini (1987) provide 
the evidence for the existence of streamwise vortices in low-speed boundary layers undergoing laminar- 
turbulent transition. Such vortices have been observed in the buffer region of the turbulent boundary layers 
by Blackwelder and Eckelmann (1979) for instance. The review of Willmarth (1975) cite some earlier 
observations of these phenomena. In free shear flows, the experimental studies (e.g., Lasheras and Choi 
(1988) and the references therein), theoretical studies (e.g., Lin and Corcos (1984) and the references cited 
therein), and the direct numerical simulations (e.g., Metcalf, et al. (1987) and Ashurst and Meiburg (1988)) 
have provided some insight into the origin and evolution of these streamwise vortices. 

The evolution of a shear flow in the presence of streamwise vortices, or conversely, the dynamics of 
vortices stretched by shear flow is of fundamental interest in fluid mechanics. A physically and geometri- 
cally simple model of such interactions between spanwise and streamwise vorticity appears to have been 
first proposed by Pearson and Abernathy (1984). Their model assumes an incompressible mean flow with 
uniform shear which is then disturbed by the introduction of an Oseen vortex with its axis in the streamwise 
direction. Corcos (1988) has studied this problem by a systematic asymptotic approach with a focus and 
range different from that of Pearson and Abernathy. He showed that the streamwise velocity component 
underwent severe distortions in a region of thickness v 1/3 as a result of the streamwise vortex. These distor- 
tions were damped out in the viscous core region, which is of size v 1/2 . This model, however, is restricted 
to the analysis of streamwise structures of turbulent non-reacting shear flows. We shall now focus our 
attention to reacting flows. 

An important problem in the general area of non-premixed combustion is the study of chemical reac- 
tions in turbulent reacting flows. However, such flows are difficult to analyze and few theoretical studies 
have been conducted. One such theoretical investigation, which forms the genesis of our model and is to be 
described in more detail later, is that of Marble (1985), who developed a simple model problem which 
locally describes an established thin flame that is being wrapped up by a small scale eddy. However, there 
are certain important aspects, such as ignition and flame structure, which were not included within the 
framework of the Marble problem. 

We have completed a systematic analysis of the ignition process and the subsequent flame structure of 
diffusion flames in vortex dominated flows using asymptotics and direct numerical simulation in the absence 
of a shear flow (Macaraeg, Jackson and Hussaini, 1992). We considered a constant-density, one-step, 
irreversible Arrhenius reaction between initially unmixed species of fuel and oxidizer occupying adjacent 
half-planes, which were then allowed to mix and react in the presence of a single point vortex. Using large 
activation energy asymptotics, particular attention was paid to the ignition regime as a function of the vortex 
Reynolds number. For near equal initial temperatures of the fuel and oxidizer and expanding the tempera- 
ture about its inert value, it was shown that the overall ignition time decreases slightly and then asymptotes 
to a finite value as the vortex Reynolds number is increased from zero. For finite values of the activation 
energy, a direct numerical solution was used to study post-ignition events. First, a hot spot develops within 
the viscous core of the vortex (ignition) and progresses rapidly into a circular flame which grows in time 
according to a similarity rule. Next, a tribrachial flame configuration evolves at infinite distance to the right 
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and left of the origin. The diffusion flames emanating from this configuration then move toward and finally 
merge with the flame in the reacted core region. We have recently extended the ignition study to flows 
with counter rotating vortices (Kozusko, Macaraeg and Jackson, 1993). 

The above mentioned problem of a flame interacting with a single point vortex is supposed to model 
certain fundamental mechanisms of turbulent diffusion flames which form the basis of many propulsion dev- 
ices. Since acoustic waves are always present in these combustion systems, the role of acoustics in such 
flame/vortex interactions was investigated (Jackson, Macaraeg and Hussaini, 1993). The acoustic field was 
considered to be a spatially uniform but time dependent pressure wave of small amplitude. This pressure 
wave was then taken to have the sinusoidal form P it) = P A sin (co t + <&), where P A is the amplitude, co is 
the frequency and O is the phase, and was regarded as a single Fourier mode of a more general pressure 
wave. It was concluded that for low frequencies, ignition is accelerated or significantly delayed depending 
on the phase, while for moderate to high frequencies, ignition is always enhanced independent of phase. 
Finally, an approximation for the ignition time was found by considering an equivalent homogeneous prob- 
lem in which spatial variations are ignored. This leads to the following implicit formula for the "homogene- 
ous" ignition time , given by 

if* 

J e p ^dx = lS, 

0 

where t£ is the ignition time in the absence of a pressure disturbance. Thus, this implicit formula relates the 
ignition time when pressure disturbances are absent to the homogeneous ignition time when pressure distur- 
bances are present. It was found that this homogeneous ignition time predicts the ignition trends very well, 
and therefore was proposed as a good engineering approximation for obtaining the ignition times for the 
more complicated flow fields. 

The above mentioned ignition studies have been confined to vortices evolving in an otherwise con- 
stant flow. However, shear is an important attribute of turbulent flows. Thus, the purpose of the present 
study is to extend the above analysis so as to investigate the effects of the interaction between a streamwise 
vortex and spanwise vorticity induced by the shear on ignition. We assume a simplified kinetic model con- 
sisting of a one-step, irreversible Arrhenius reaction between unmixed species occupying adjacent half- 
planes which are then allowed to mix and react in the presence of a streamwise vortex embedded in a shear 
flow. We remark here that any realistic modeling of complex kinetics will necessitate a full numerical solu- 
tion. Then questions arise as to the reaction rates and their dependence on temperature and the relative 
importance of various reactions which entail thus rather large uncertainties. In a mathematical treatment of 
combustion, the kinetic model has to be necessarily simple. The one-step irreversible Arrhenius model has 
been extensively used with significant success in the study of low-speed combustion (e.g., Buckmaster and 
Ludford, 1982). This model appears to cover the essential physics of the problem. In the last decade or so, 
the asymptotic studies of combustion based on this model have significantly enhanced our understanding of 
ignition, of flame stability and of diffusion flame structure in low subsonic flows. There is as yet no reason 
to believe that this model will not play an equally significant role in enhancing our understanding of super- 
sonic combustion. In any event, this idealization makes the problem amenable to asymptotic analysis and 
thus provides a semi-analytical solution. The results can then certainly verify and in turn be verified by full 
numerical simulations. An idealized model is now employed to describe the interaction process: the 
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streamwise vortex, assumed to be infinite in extent, interacts in time with the streamwise velocity com- 
ponent associated with the spanwise vorticity. Initially, the streamwise velocity is assumed to be a linear 
shear profile and the streamwise vortex resides along the fuel-oxidant interface. The reactants are then 
allowed to mix and react. Emphasis is placed on ignition time and structure. In Section 2 the basic equa- 
tions and assumptions are presented. Section 3 contains the ignition analysis and results, and the conclusion 
is given in Section 4. 

2. PROBLEM FORMULATION 

In this section the problem for the time evolution of initially unmixed species occupying adjacent 
half-planes which are then allowed to mix and react in the presence of a vortex embedded in a shear flow 
(Figure 1) is formulated. We assume that the flow is independent of the downstream direction x for all 
time. Although it is not strictly true for any real flow, it is consistent with experimental observations that 
streamwise vortical structures in turbulent flows are essentially aligned with the flow and are of great 
streamwise extent Since we are interested in flows with vortical motion, the basic equations are 
transformed into cylindrical coordinates (y ,z) -> (r,0) by the transformations y = r cosG and z = r sinO. 
If we define g to be the radial velocity and h to be the axial velocity, i.e., v = g cosG - h sin 6 and 
w = g sin9 + h cos 8, then the governing non-dimensional equations in cylindrical coordinates, neglecting 
all dependency in x , reduce to 


pT =P , 


Pr + (p#)r + “(p h)o + = °* 


(2.1a) 

(2.1b) 


p[M< +8 U r + j Uq] = ^ VV 


(2.1c) 


r A h 2 , p r 1 

PU. + S&- + -S9-— ] + ^2 = ^r 




r r 


jh + jK r 


(2. Id) 


(,lhl+g hr + jh e + -^-]+-y^--^ 


V 2 A-4 + 4 Se+T-^9 

r 2 r 2 3 r 


CUe) 


P iT t 


\+ g T r + ^T']-^lP l+g P r + ^P e ] = j^V 2 T + !X- 1 ^<b + VCl, (2.1f) 


Re 


p[Fj.,+gFj.r + jFj,,]= -O. J = ‘.2. 


(2.1g) 


Q = t Da pFiF 2 e 


-ZelT 


(2.1h) 
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where V 2 is the two-dimensional Laplacian operator in cylindrical coordinates, is the viscous dissipation 
term, and K = g r + g / r + h e / r . Here, (u,g,h) are the velocity components in the (x , r, 0) directions, 
respectively; p is the density; T is the temperature; P is the pressure; and F { and F 2 the mass fractions of 
the fuel and oxidizer, respectively. The chemical model is assumed to be a one-step, irreversible Arrhenius 
reaction. The nondimensional parameters appearing above are the Schmidt number Sc =v/D assumed 
equal for both species with D the species diffusion coefficient and v the kinematic viscosity; the Prandtl 
number Pr\ the Reynolds number R t = L v VJv where L v and V v are characteristic length and velocity 
scales of the vortex; the Zeldovich number Ze =E /(R° rj with E the dimensional activation energy and 
r° the universal gas constant; the Damkohler number Da = t d /t c defined as the ratio of the characteristic 
diffusion time scale t d to the characteristic chemical time scale t c ; p the heat release per unit mass of F y 
the ratio of specific heats; the Mach number M = V v / a „ where is the characteristic speed of sound; and 
finally x = t v / t d is the ratio of the time scale of the vortex t v to that of diffusion. The density, temperature, 
pressure and mass fractions are nondimensionalized by their initial values p M , T rM> / > eo = p co /?°r oo and F l oo , 
respectively. Lengths are referred to L v , a characteristic lengthscale of the vortex, such as the core size, and 
velocities are referred to V v , the reference tangential velocity. The actual choice for these scales is given in 
the next section. The relevant time scale is therefore t v = L v /V v . The system (2.1) must be solved subject 
to appropriate initial and boundary conditions chosen to represent a streamwise Oseen vortex evolving in a 
shear flow. 

3. IGNITION 

At time t - 0, the reaction rate is exactly zero owing to the product F 1 F 2 = 0. For t > 0, the fuel and 
oxidizer begin to mix by diffusion, as well as by convection due to the presence of the vortex, and the reac- 
tion rate is no longer zero. For small time, it can thus reasonably be assumed that the effect of the reaction 
on the overall flow field is small. The solution for small time in which the reaction rate term is neglected is 
known as the inert or chemically frozen solution, and will be denoted by the superscript I. In what follows, 
we only consider the case of near equal initial temperatures and small Mach numbers, and take the asymp- 
totic limit of large Zeldovich number (for the simplified hydrogen-oxygen system 2H 2 + 0 2 — > 2H 2 0 in 
the presence of an inert N 2 , the dimensional activation energy E lies in the range 0(7-40 kcal/mole) yield- 
ing Ze -5-20 for in the range 300 - 2000° K; e.g., Glassman, 1977; also, Drummond and Mukunda, 
1988, in their numerical work on reacting mixing layers uses 7200 cal/mole). For large activation energies 
(and hence large Zeldovich numbers), this assumption leads to the inert solutions T 1 = p 1 = 1 + O {Ze~ l ), 
provided that the Mach number is no larger than Ze~ m so that the viscous dissipation term can be neglected 
to leading order. For larger Mach numbers the temperature field can no longer be considered uniform to 
leading order since heating must take place by means of the dissipation term. The 0(Ze~ { ) term is included 
to allow for small initial temperature differences across the interface. Thus, the leading-order inert solution 
corresponds to the constant density approximation prevalent in the combustion literature. As time increases, 
more of the combustible mixes until, at some finite time, a thermal explosion occurs characterized by 
significant departure from the inert. To analyze the ignition process, we determine the effect of the growing 
reaction rate by expanding about the inert solution as 

T - 1 + Ze" 1 T { + O (Z<T 2 ), p = 1 + Ze~ 1 p l + 0 (Z<T 2 ), F } = F) + O ( Ze ~ x ), (3.1a) 
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(u,g ,h) = (u 1 ,g l ,h ! ) + 0(Ze- 1 ), P =1 + yM 2 [P l +0(Ze x )], (3.1b) 

and take the asymptotic limit Ze — > =*>. The viscous dissipation term in (2. If) in this limit is of the order 
M 2 Ze , and for a proper balance of the energy equation we set 

M=M 0 HZe. (3.2) 


Thus, solutions are also sought in limit of small Mach numbers. A comment about the size of the Mach 
number and the relative importance of the viscous dissipation term in the energy equation is in order here. 
For typical combustion systems the Mach numbers are of the order O(10 -3 ), and so the role of viscous heat- 
ing is negligible. Here, however, we have in mind propulsion devices in high speed flows. As such, the 
Mach number within the mixing region can be relatively large. For example, in their study of compressible 
viscous vortices, Colonius et al. (1991) report Mach numbers in the range of 0.0625 to 0.67 (the definition 
of their Mach number is the same as used here). In our work, we have taken M 0 = 'IZe M with M 0 lying 
in the range (0.3-0.7), corresponding to Mach numbers of 0.01 to 0.22, respectively, for a typical value of 
Ze=10. These values are well within physical ranges. Thus, in these flow situations, viscous heating can 
indeed play a major role in the ignition dynamics. 

The leading-order equations are given by 

Pl + T x = yM 2 P * , (3.3a) 

g! + jhl + ^r = 0, (3.3b) 


u! + g ! 4 + -y- u l e = ^ W, 


(3.3c) 


g! + g 1 g l r 


+ hL A -h!L + P i = ±. 

r 50 r r Re 


vy -y - 


sL 2^ h i 

T^e 


(3.3d) 


*/ + ,'A' + 4*i + i T L + 7 , ’ t "to 


vy - 4 + 


r 2 r 2 


(3.3e) 


T Ut +g l 


1 T Ur + -T U9 - (Y- 1 )M 2 P! = V 2 ^ + (Y~ 1 ) ^ & + xF} Fi e 


(3.3f) 


FL+jF! r + 


J't 


hL 

r 


F/.e 


Re Sc 


V 2 Fj , y = l,2. 


(3.3g) 






(3.3h) 
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where we have chosen the Damkohler number Da to be 


Da 


,Ze 


zfiZe 


(3.4) 


This particular choice of the Damkohler number ensures that a distinguished limit exists, in that the reaction 
rate term is of the same order in a Zeldovich number expansion as the time derivative terms (see, e.g., 
Buckmaster and Ludford, 1982). This implies that we are looking at the fast chemistry limit and thus selects 
our chemical time scale. 

We first see that the continuity and momentum equations (3.3b-e) are decoupled from the rest of the 
system and thus, will completely determine u 1 , g 1 , h 1 and P 1 . At t = 0, we situate a point vortex at the ori- 
gin and allow it to diffuse under the action of viscosity. Then, assuming the flow in the (r,0) plane to be 
axisymmetric, the induced non-dimensional flow field is given by (Lamb, 1932) 


g l = 0 , h 1 = —■ [1 - exp(-Rer 2 /4 0L 


(3.5) 


P l (r ,t) = -J Vfe iil dp = j ~~ (1 - 2e~*' 4 + e~^ 12 ) + Ei (t\ 2 /4) - £/(n 2 /2) j , (3.6) 

where r \ = r VRe/VF. The solution (3.5)-(3.6) is the solution for the incompressible Oseen vortex. The 
velocity scale V v has been chosen so that the corresponding potential vortex at L v (i.e., at r = 1) is unity, 
resulting in the relationship 


27 iL v Vv 


= 1 , 


where F is the circulation parameter of the vortex. This particular choice shows that 

r 


Re = 


L v V v 


2nv 


= R« 


where R a is the usual vortex Reynolds number. The stream wise velocity component u 1 is found from the 
equation 


h l 




(3.7) 


which is to be solved subject to the initial and boundary conditions, 

u I =k i + k 2 z at t > 0, r 0 < 0 < 2jc, and t = 0, r > 0, 0 < 0 < 2tc, (3.8) 

where ki is the nondimensional velocity at the centerline and k 2 is the local nondimensional shear, defined 
as the ratio of the time scale of the vortex to that of the shear. In this study we shall assume that k 2 ~ O (1) 
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so that the two time scales are comparable. At this order the streamwise vortex influences the streamwise 
velocity component and hence the associated spanwise vorticity, but not vice versa. 

We note here that equation (3.7) has previously been considered by several authors for non-reacting 
flows. Pearson and Abernathy (1984), in the context of streamwise vortical structures found in turbulent 
boundary layers, constructed solutions which were valid for all Reynolds numbers. Corcos (1988), in the 
context of large scale structures found in mixing layers, has recently given an asymptotic solution in the 
limit of large Reynolds number. Because these authors used different notations and because these solutions 
are needed later for understanding certain aspects of ignition, we give below a summary of these solutions 
for completeness. 

To examine the time evolution of the streamwise Oseen vortex in a shear flow, we define 

« J (r,9,0 = ki + k 2 [r sin9 + u(r,Q,t)] , (3.9) 

where m is the disturbed velocity component owing to the presence of the vortex. Upon substituting (3.9) 
into (3.7), u is found to satisfy the convection-diffusion equation 

h, + \ [1 - exp(-Rer 2 /4f)] (r cos0 + n e ) = V 2 /?, (3.10) 

K.C 

subject to homogeneous initial and boundary conditions. This equation was previously considered by Pear- 
son and Abernathy (1984) (see their equation (14)), in the context of streamwise vortical structures found in 
turbulent boundary layers. They employed the similarity type variable r 2 /4t and defined (in our notation) 
u = Real(r F(r)e i6 ). This choice led to an equation for the complex function F(r) which was solved 
numerically after special treatment of the regular singular point at the origin. A different approach for solv- 
ing (3.10), which is perhaps more instructive, is to first transform using the similarity variable 
T| = r VRe/V7 and then introducing the definition 

u = —L- [u e (T|)cos0 + u 0 (n)sin0]; (3.11) 

\Re 

this yields the coupled second-order differential equations, 

— (m c — 1) ) h — 2 [1 — exp(— T| 2 / 4)J (u 0 + Tj ) = u e + — u e — 2 u e > (3.12a) 

l T| T 1 M 

-J- ( u 0 - 11 u ' 0 ) - ^ [1 - exp( — rj 2 / 4)] u e = u" + ^ u' 0 — ~ u <> > ( 3 - 12b ) 

2 Tj M «1 

subject to the boundary conditions, 

u e (0) = u 0 (0) = u e H = u 0 (<~) = 0. (3.13) 
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Primes denote ordinary differentiation with respect to the similarity variable T|, u e is the amplitude of the 
even mode, and u 0 is the amplitude of the odd mode. The system (3.12) depends only on the Reynolds 
number Re. 

The coupled system (3.12) for the odd and even components subject to (3.13) represents the influence 
of the Oseen vortex on the streamwise velocity and is solved numerically using a 2nd-order finite difference 
scheme with an appropriate stretching in T|. The solution to this system for a range of Reynolds numbers is 
given in Figure 2. As the Reynolds number increases, the number of oscillations in each of the components 
u e and u 0 increases while the region of oscillation moves away from the origin, thus establishing a core 
region. The solution difference u d = T| sin 0 + u e cos 6 + u a sin 0, where u d - ^lR&(u r - ki)l(k 2 ^t ), is 
plotted in Figure 3 for 0 = ±n/2 and for various values of Re. From this figure the extent of the core 
region is seen to grow as Re increases, and this core, once established, grows like FF . 

To investigate the solution described above for large times, let t -> <», so that rj -» 0 for fixed r . 
Hence from (3.11) we have 

<t 

u - lim -t= [« c (t|)cos0 + w o (n)sin0] 
r vRe 

= lim — [M e (n)cos0 + K<,(n)sin0] 

T| — >0 TJ 

- lim + K„0l)z] 

Ti — >0 T| 

= aiy+^iz, (3.14) 


where from (3.12) and (3.13) it is easy to show that the asymptotic behaviour as t| — > 0 is 

u e - a\X\ + 0(il 3 ), u 0 =bi\\ + 0(X\ 3 ), (3.15) 

where a\ and b\ are found numerically and are functions of the Reynolds number Re. Substituting (3.14) 
into (3.9) yields 


u 1 = ki + k 2 [a 1 y + (1 + &i)z ]. (3.16) 

Comparing (3.16) to the initial condition (3.8), the effect of the vortex is seen to be the establishment of a 
new equilibrium shear profile within the viscous core region. In addition, as Re — > <*> the coefficients 
a i -> 0, bi -» - 1, yielding u 1 -> i.e., the streamwise velocity component becomes constant within the 

core at large vortex Reynolds numbers. 

The streamwise Oseen vortex induces additional spanwise vorticity. This is seen by considering the 
total vorticity defined by © = e x co + m , where 0 ) is the streamwise vorticity component and © is associated 
with the spanwise vorticity 


© = 


„ du 

e > aT 



(3.17) 
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The magnitude of the spanwise vorticity, defined by G = | © \ , is a function of Re and, at leading order, 
has the value 


Go - | &2 I 



+ (1 + b\) 2 



(3.18) 


at the center of the core. Figure 4 is a plot of G 0 versus Re for k 2 = 1. Note that as Re -> G 0 —> 0 
which implies that there are no gradients within the center of the core and thus the core is completely esta- 
blished. From this figure, the core is 95% established at Re = 50 which is in agreement with Figure 3 where 
a well defined core is clearly visible at Re = 100. 

Even though the solution can be determined numerically for any Reynolds number, it is constructive 
to consider the limit as Re -> °°. The leading-order solution in this limit is found by a direct asymptotic 
investigation of (3.7) subject to the initial and boundary conditions (3.8). In the limit Re -> an outer 
inviscid layer exists and is governed by the inviscid version of (3.7), 

«,°+4«e= 0, (3.19) 

r L 

where u° denotes the streamwise velocity u l in this outer layer. This equation is valid for r - (7(1) and 
t < 0(Re ), so that the effect of a potential line vortex or an Oseen vortex on the streamwise velocity is the 
same. The solution found by the method of characteristics is 

u° = ki + k 2 r sin 0 O , (3.20) 


where 9 q is the characteristic 


e 0 = e - 4- < 3 - 21 > 

r 

As r 0, the solution (3.20) is singular indicating that the viscous terms in (3.7) must be retained. 
We introduce the similarity type variable 



( 3 . 22 ) 


where Re -1 <<8<<1 and 8 is a nondimensional parameter that will be chosen in the course of the 
analysis. Defining u M to be the dependent variable u l in the main layer, and letting 


u M = k i + £2^5* n(s,0). 


(3.23) 


we see that u satisfies the convection-diffusion equation 

Vi(u -h s ) + -tt [1 ~ exp(-ReSs 2 /4)] m 0 = -r—— V 2 n. (3.24) 

os oKe 
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To facilitate matching with the outer inviscid solution (3.20), we change to a new (non-orthogonal) coordi- 
nate system in which the characteristic 9 0 is one of the coordinates. Thus, the above equation in terms of 5 
and 0 O becomes 


Vi(u - 5 u s ) - — -exp(-ReSj 2 /4)u e 

8.r 0 


a 2 a 


5 Re 


a s 8 s 3 aa 0 


i 

+ — 
s 


d_ 2 a 

ds bs 3 a©o 


i 2 39o 2 | (M> ' 


(3.25) 


Balancing the viscous terms with the convection terms necessitates the choice 


i = _i_ A. 


(3.26) 


or, in terms of 8, 


8 = (4/Re) m s 0(v 1/3 ), (3.27) 

which is in agreement with Corcos (1988). This layer is thicker than the inner diffusion layer of the vortex 
which is of the order v 1/2 . For this scaling, exp(-Re8s 2 /4) is exponentially small in the distinguished 
limit Re >> 1, and so again the effect of a potential line vortex or an Oseen vortex is the same. In order to 
account for the diffusion of the vortex, an inner layer must be constructed. This is not necessary, as will be 
shown below, since u M is not singular at the origin. The leading order equation in a 8-expansion is given 
by 


V2(u - su s )= “ «0 o e o , 


(3.28) 


which is solved by separation of variables and matched with the outer solution (3.20) to get 


u = s exp 


-1 

3s 6 


sin0r 


(3.29) 


-1 

3j 6 


sin 0r 


The solution in the main layer is now given by 

u M = ki + k 2 s *f§7 exp 
The composite solution, valid in both the outer and main layers, is given by 

u c = ki + k 2 r exp 


(3.30) 


’-(S/) 3 ’ 

sin 

* ' 

0- -y 

3r 6 


r 2 


(3.31) 
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To compare the asymptotic solution u c given in (3.31) to the exact solution u l given in (3.9), a plot of the 
solution difference u d versus T| for 9 = ±%/2 with Re = 10,000 is shown in Figure 5, where u d and T| are 
defined as before. From this figure we see that the asymptotic solution is in good quantitative agreement 
with the exact solution shown in Figure 3c. 

The magnitude of the spanwise vorticity at 0( 1) obtained by substituting (3.31) into the definition of 
G and then letting 5 -» 0, yields 

cos(0-— -r)i. (3.32) 

or 


G 0 = 


2[ * 2 | 

&s 2 


exp 


-1 

3s 6 


Figure 6 is a plot of G 0 versus s for 0 = rc/2, k 2 = 1, and 8 = 1/5. The limit Re -» 0 corresponds to 
Go — > 0, and so the large spikes have been created by the vortex. 

With g ! , h 1 and u 1 now known, the following equations can be solved for T\ and F- 


i ,t 


h 1 

1 1,9 
r 


R QPr 


V 2 rj + (t- l)M 0 2 [P! +~^]+ Fi e T \ 


Re 


(3.33) 


17/ 1 Jl pi 

F j,t + - F j,Q 


Re Sc 


V 2 Fj, j = 1,2. 


(3.34) 


Note that the right hand side of the equation for the temperature is the sum of three terms: the diffusion 
term, the term due to compressible effects, and the reaction term, respectively. The appropriate boundary 
and initial conditions are given by 

T 1 = 0, F[ = 1, F l 2 = 0 at t = 0, r > 0, 0 < 0 < n, and t > 0, r -> 0 < 0 < n, (3.35a) 


T x = p r , F[ = 0, F l 2 = (J) -1 at t = 0, r > 0, n < 0 < 2n , and t > 0, r -» % < 0 < 2 it, (3.35b) 

where is the equivalence ratio defined as the ratio of the initial mass fraction of the fuel F lt „ to the initial 
mass fraction of the oxidizer F 2i _~, and is the parameter which allows for small initial temperature 
differences across the interface. If <j) = 1, the mixture is said to be stoichiometric; if <|> > 1 it is fuel rich; 
and if $ < 1, it is fuel lean. Also, if is less than zero, the oxidizer (species 2) is relatively cold compared 
to the fuel (species 1); and if p T is greater than zero, it is relatively hot. As t increases, the solution for T { 
becomes unbounded at some finite time ( t ig ) and location (y ig , z ig ). This characterizes the ignition regime. 
Finally, with T v determined, the density perturbation p! can be found from (3.3a). 

The case M 0 = 0 was discussed in the Introduction. Thus, the goal is to examine the influence of the 
compressible term in (3.33) on ignition. However, it is instructive to first investigate its influence on the 
temperature field in the absence of chemistry. We do this in the next section. The following section then 
contains numerical results when chemistry is included. 

3.1. Nonreacting Case 
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In this section we investigate the influence of the compressible term in (3.33) on the temperature field 
in the absence of chemistry. Since an exact solution of (3.33) with t = 0 is not possible, we shall examine 
the equation in the limit Re >> 1. In this limit an outer inviscid layer exists and is governed by the inviscid 
version of (3.33) for the outer variable T ? . Transforming to the (r,0 o ,f) characteristic coordinate system, 
where 0 O is given in (3.21), the temperature equation (3.33) becomes to leading order 


T 


o _ 
i,< - 


1 

R ePr 


ff_a_ it a 

r 3 38 0 



d 2 1_ d 

dr r 3 00 q 


i a 2 

r 2 30O 2 


(Tf) 


+ 


7- 1 
Re 


M„ 



U r 



+ ~T &8 n 


+ 0 (Re ~ 2 ) , 


(3.36) 


where u = r sin 0 O and we have made use of the leading order outer solutions h 1 = Hr + est and 
p 1 = - 1/(2 r 2 ) + est, where est means exponentially small terms in Re, valid for t « O (Re). The lead- 
ing order solution in a large Reynolds number expansion can be shown to be 


T? = ^M 3 



t 2 4t 3 

t + — sin(20 o ) + — j 
r 3 r 


cos 1 0 O 


+ O {Re- 2 ) 


(3.37) 


Note that the terms in the brackets are associated with the streamwise vortex and with the spanwise vorti- 
city. In the absence of the spanwise vorticity (i.e., = 0), the first term is positive showing that the effect 

of compressibility is to heat the outer region of the viscous core. With rf known, the density in the outer 
layer can be determined from (3.3a). 

As r -» 0, the solution (3.37) becomes singular so that the viscous terms in (3.33) must be retained. 
Transforming to the (s,0 o ,O coordinate system, where s is defined in (3.22) and 0 O = 0 - 1 /Sj 2 , the tem- 
perature equation (3.33) in the main layer becomes 



-V 'izr yaf) 


(3.38) 


where, to leading order in 6, 


u = s exp 



sin0 o + 0(8). 


(3.39) 


To determine the proper expansion in the main layer and the appropriate matching condition as s -> we 
first expand the outer solution (3.37) as r -» 0, yielding 


T\ ~ (Y - 1 )M 2 


kit 2 


w b 2 kit b k% t 

1 + — - — cos 2 0o + sin (2 0 O ) + — — — 

4 s 4 


; 3 fc 2 


t s 


+ 


(3.40) 
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We now expand Tf as 


rf = 8/i + 8 2 / 2 + 


(3.41) 


as 8 0. Substituting the above expansion into (3.38) gives the following problem at leading order: 

-2 


*2 


f, - — f , = — !— /\ + (Y- 1)M 2 J + — exp 

7l ‘ 2f 71 ' Vo U 1 fV j 4 


3s f 


cos 0 O f . 


(3.42) 


The appropriate solution that matches with the outer solution (3.40) can be verified to be 


CY-l) Mj Tr i 
fl ~ ti* 4 


Mo ki t f n(5 ,0o)> 


fn = * 2 i 1 - exp 


3s t 


+ 


Pr 

(2 -Pr) 


exp 


-2 

3s 6 


- exp 


-4 

3Pr s 6 


cos(20 o ) > , 


(3.43a) 

(3.43b) 


which decays as s -> *». The first term of /j is associated with the streamwise vortex while f n is associ- 
ated with the spanwise vorticity. As mentioned above, the leading order term f x is positive showing that 
the effect of compressibility is to heat the outer region of the viscous core. Since we are primarily 
interested in the influence of the spanwise vorticity, we plot in Figure 7 a graph of / x x as a function of s 
for Pr = 1, 0 = nil, and (a) 8 = 1/5 and (b) 8 = 1/8. Note that the graph decays exponentially as s -» «>, 
quadratically as s -> 0, and is positive and oscillates for s ~ 0(1). Also, the number of oscillations 
increases as 8 decreases but stays within some temperature envelope. Thus the effect of the streamwise vor- 
tex is to introduce positive temperature fluctuations in the main layer which are damped out in the inner 
viscous layer. This result will have important implications in the ignition regime when chemistry is 
included, and will be brought out in the next section. The density can now be determined from (3.3a), and 
is given by 


pi“— 


■jMj 
2 s 2 t b 


-8/1 + 


(3.44) 


As s — » 0, rf and pf* become singular suggesting that an inner viscous layer is present. To deter- 
mine the proper scalings, we see from the expansion (3.44) for the density that the second order term is of 
the same order as the leading order term when s 0, yielding the balance s = 0(Re -1/3 ). This suggests 
the new variable for the inner layer 



(3.45) 


which is the size of the viscous core. Expanding the temperature and density in the main layer in terms of 
the inner variable and taking the limit Re -> 00 , the proper expansion for the inner variables are 

T[ - Ref, pi = Re p, (3.46) 
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where the superscript indicates variables in the inner layer. The leading order equations in a Reynolds 
number expansion are given by 


p + T = yM?P 


(3.47a) 


f, - (y- 1 )M„ 2 F 0 , = + Y ) + (V- - j- )2 


(3.47b) 


where 


h 0 = Rz m h 0 ttj), P o - ReP 0 (tO; (3-48) 

fiQ and P 0 can be found from (3.5) and (3.6). This system is axisymmetric and is to be solved subject to 
appropriate matching conditions as ^ We do not solve this system here, but rather note that the vari- 
ables are bounded at the origin. Furthermore, this system is equivalent to that considered by Colonius et al. 
(1991) if we rewrite their system in our notation, and so we refer the interested reader to that paper for 
further details. 

3.2. Reacting Case 

When chemical effects are included, the system (3.33)-(3.34) with boundary conditions (3.35) 
describes the ignition regime which must be solved numerically. The numerical scheme utilized here is 
described in Macaraeg et al. (1992). A small nonzero initial time is assumed in order to avoid nonphysical 
initial transients which may lead to early ignition. For the cases which follow, k\ is assumed to be zero, so 
that the imposed shear is symmetric about the origin. Emphasis will be placed on the effects of the shear 
strength parameter k 2 , so the other parameters will be held fixed. Here, we set = 0 and (j> = Pr = 1. 

We begin the numerical computations by first setting t = 0 in (3.33). This allows a qualitative com- 
parison between the numerically generated solution and the asymptotic solution derived in the previous sec- 
tion. To this end, recall that the asymptotic solution indicates that the introduction of streamwise shear 
gives rise to oscillations in the velocity field as well as the temperature field via the dissipation function. 
Shown in Figure 8 are contour plots and cross-section cuts of T i and TU for Re = 50, k 2 = 1 and M 0 = 0.7. 
The quantity TU is computed by subtracting the numerical solution corresponding to M 0 =0 from TV Thus, 
TU is proportional to /n of the previous section and highlights the underlying oscillating temperature field. 
The cross-section plots of T x and TU versus z , taken through the vortex center, show that the overall tem- 
perature perturbation has a single temperature peak in the center, whereas the temperature due to spanwise 
vorticity has multiple peaks, even though their magnitude is small at this low shear strength. If the strength 
of the shear is increased, the magnitude of these oscillations increases by k 2 , as can be seen in Figure 9 
which displays plots of TU for three values of the shear strength parameter k 2 . This trend was predicted 
earlier by the outer asymptotic solution obtained in (3.43), which indicates a proportionality with k 2 . Also 
given are the corresponding cross-sections of the streamwise velocity component. Note that the magnitude 
and number of oscillations increases with increasing shear, again consistent with the asymptotic solution 
derived in the previous section. Figure 10 depicts contours and cross-sections of T t for the same conditions 
of Figure 9. Note that the temperature perturbation, TV becomes dominated by the spanwise vorticity for 
high values of shear. 
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The multiple temperature peaks appearing for high shear values can have a dramatic effect on ignition 
when chemistry is added. To illustrate this we show in Figure 11 contours and cross-section plots of the 
temperature perturbation for Re = 10, M 0 = 0.316 and two values of the shear strength parameter ft* For 
*, = i, there is a single hot spot located in the vortex center, while for high shear (ft 2 - 20) there are two 
ignition points, located in the main layer of the viscous core. The evolution in time of this double ignition 
point is shown in Figure 12. Note that the temperature perturbation is strongly sheared, which gives rise to 

multiple ignition points at a later time. 

The influence of the streamwise shear on the ignition time is shown in Figure 13, where the maximum 
temperature perturbation over the entire spatial domain is plotted against time. We note here that for weak 
shear (ft 2 = 1), the ignition time is comparable to that obtained by Linan and Crespo (1976) m which 
diffusion is the only mechanism present for mixing. Also, the ignition location is at the center of the vortex, 
consistent with previous results (Macaraeg et al„ 1992). As the strength of the shear increases, the ignition 
time decreases indicating that mixing is enhanced owing to the presence of streamwise shear. In addition, 
the ignition location no longer resides in the viscous core, but now multiple igmuon points occur and are 
located in the outer region of the core. This phenomena is due to the large temperature perturbations 
induced by the interaction process between the streamwise and spanwise vorticity described in the previous 

section. 

4. CONCLUSIONS 

The problem of a flame interacting with a vortex is assumed to model certain fundamental mechan- 
isms of turbulent diffusion flames which form the basis of many propulsion devices. Streamwise shear, 
which can be present due to large scale structures within the flow field or can be introduced as a mixing 
enhancement technique, is always present in these combustion systems, so that it is important to study the 
role of shear in such flame/vortex interactions. An unexplored fundamental problem of turbulent reactmg 
flows is thus the evolution of a shear rfow in the presence of streamwise vortices, or conversely, the dynam- 
ics of vortices stretched by shear flow. This problem is examined by means of asymptouc analysis and 
numerical simulation in the limit of small Mach number. An idealized model is employed to describe toe 
interaction process which, among other assumptions, neglects spatial variations in toe streamwise direction. 
The model then consists of a one-step, irreversible Arrhenuis reaction between initially unmixed species 
occupying adjacent half-planes which are then allowed to mix and react in toe presence of a streamw.se 

vortex embedded in a shear flow. 

It is shown that the effect of shear is to give rise to oscillations in the velocity field, which produces 
localized regions of high temperature gradients via viscous heating. As toe strength of toe shear is increased, 
these high temperature gradients give rise to multiple ignition points while substantially decreasing toe 
overall ignition times. This results in the enhancement of mixing brought about by the linear shear associ- 

ated with the spanwise vorticity. 
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2 . Plot of (a) u e versus 77 for various values of the Reynolds number Re. 
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5. Plot of the solution difference Ud versus rj for 0 — ±7r/2 and Re = 10,000. 
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. Plot of /n as a function of s with Pr = 1,0 = tt/2 and for (a)6 — 1/5 and 8 — 1/8. 






9. Vertical cuts at the vortex center of temperature perturbation due to the spanwise vorticity 
(TU) and streamwise velocity (£/), for varying shear strengths. (No chemistry) Re = 50, 




R 0 = 10 M 0 = .316 




12. Contours of the temperature perturbation (Ti) evolution from early time to ignition. 
Re = 10, M 0 = 0.316, k 2 = 20. 
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